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A triple system with a binary on the same line of sight 
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Abstract. UX Ari belongs to the class of very active RS CVn stars and has recently been the target of surface 
(Doppler) imaging. Although this technique needs a quite accurate determination of the orbit (in order to have 
the correct period for phasing and the correct Doppler shift correction of the line profiles) we found only one, 
quite old orbit solution, which has subsequently been used by everyone. 

We used published radial velocities (RVs), supplemented by a large number (124) of our own recent, high-accuracy 
RVs of both the primary (KG IV) and the secondary (G5 V) to improve the orbit of UX Ari. In addition to the 
improved set of parameters, we found that the 7 velocity of the system is systematically changing over time. It 
seems that UX Ari is a triple system. Actually, a third star is weakly present in the spectrum. While its RV is 
also changing, it is not a member of the system, but happens to be on the same line of sight. 
Finally, conclusions about the physical parameters of the objects from the orbits are presented. 
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1. Introduction 

UX Ari (HD 21242; KOIV + G5V) is a short-period 
{P « 6^4), double- lined spectroscopic binary. It belongs 
to the class of RS CVn stars, i.e. at least the cool primary 
shows signs of activity. Thus, it is listed in the catalogue 
by Strassmeier et al. ( 19931 ), where more information on 
the system may be found (UX Ari=CABS 28); it was also 
put in our long term programme of surface im aging of 



1998) 



active stars (for results see e.g. Berdyugina et al. 
Surface imaging needs good orbital parameters. The spec- 
tral lines, whose distortions are followed through the ro- 
tational phases^ need to have the Doppler shift due to the 
orbital motion removed; while experience shows that small 
residual shifts do not change the main surface structures. 



Send offprint requests to: R. Duemmler 
e-mail: Rudolf .Duemmler@Oulu.Fi 

* based on observations collected at the Nordic Optical 
Telescope (NOT), European Northern Observatory, La Palma, 
Spain, and at the McDonald and Kitt Peak Observatories, USA 
** Table 2 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr ( 130.79.128.5) or via 



http:/ /cdsweb. u-strasbg.fr/ Abstract. html 



^ given the evolved state of one of the components and the 
short period, it is reasonable to assume synchronized rotation. 
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the extra noise in the data might lead to a lower quality 
of the map. Normally, the cross-correlation technique can 
be used to align spectra and remove radial velocity shifts 
without knowledge of any orbital parameters. However, 
cross-correlation relies on the assumption that the spectral 
features in the programme spectrum and the template are 
identical, which they are not in the case of active stars due 
to the distortions caused by spots. This leads to system- 
atic radial velocity errors. It is hoped that if one compiles 
a large set of radial velocities spanning a long time, the 
constraints of orbital motion and the limited lifetime of 
spots allow the determination of a good set of orbital pa- 
rameters despite the fact that individual radial velocities 
obtained in a short time span are systematically shifted. 
In order to calculate the rotational phases, a good value 
for the period needs to be known. This is especially im- 
portant if maps obtained during several seasons are to be 
compared: an incorrect period and conjunction time lead 
to increasing phase shifts which mimic motions of the sur- 
face structure which are not real and lead to incorrect 
interpretation of the long term behaviour of the surface 
structures. Thus a (re-)determination of the orbit prior to 
surface imaging is strongly recommended. 

For UX Ari, there seems to exist o nly on e orbit determi- 
nation: that by Carlos & Popper ( 1971 ). Given the age 
and the low number of measurements used, we felt it long 
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Table 1. The wavelength resolution {R = A/AA) and 
range in signal-to-noise ratios (S/N) at 6400 A for each 
dataset. 



Set 


R 


S/N 


S95 


86 000 


270-420 


S96 


71000 


280-390 


M99 


36 000 


180-390 


K99 


86 000 


90-170 


MOO 


48 000 


140-360 



Table 3. Measurements from speckle interferometry and 
from Hipparcos (ESA 1997) of the angular separation i? 
between the RS CVn system and the third star at different 
epochs (as Besselian year). 



Epoch 


^(") 


1985.8431 


0.432 


1991.25 


0.340 


1995.9237 


0.297 


1996.8658 


0.256 



Reference 



McAlister et al. (|1987| ) 
Hipparcos (ESA 



Hartkopf et al, 
Hartkopf et al. 



1997 



( 1997 ) 
(200C) 



overdue to compute a new orbit and improve the orbital 
parameters as much as possible. The time difference be- 
tween the first observation given by Carlos & Popper and 
our last observation is more than 42 years, which lets us 
hope to significantly improve the period and subsequently 
all other orbital parameters. 

Another finding makes an accurate inspection of the orbit 
particularly interesting. Lestrade et al. ( 1999 ) performed 
high-precision VLBI astrometry of, among others, UX Ari. 
They used observations at 10 epochs between July 1983 
and May 1994 to calculate the acceleration of UX Ari's 
position in the sky. The resulting acceleration was much 
larger than the perspective secular change in proper mo- 
tion and could be caused by a third body. According to 
Lestrade et al. ( 1999D , the orbital period of this system 
should be many times their 1 1-year VLBI data span. 



2. Observations and data reductions 

The new high-resolution, high signal-to-noise spectra were 
obtained for the purpose of surface imaging during five 
observing runs using three telescope-instrument combina- 
tions: 

— 5 spectra obtained in December 1995 using the high- 
resolution echelle spectrograph SOFIN (Tuominen 



1992) mounted at the Cassegrain focus of the 2.56 m 



Nordic Optical Telescope (NOT) on La Palma, Canary 
Islands, Spain (hereafter denoted S95) 

— 10 spectra obtained in November/December 1996 us- 
ing NOT/SOFIN (hereafter denoted S96) 

— 64 spectra obtained in January 1999 using the 
Sandiford Cassegrain Echelle Spectrograph mounted 
on the 2.1m Otto Struve Telescope at McDonald 
Observatory, Texas, USA (hereafter denoted M99) 

— 23 spectra obtained in February 1999 using the 
Coude CCD Spectrograph receiving light from the 
0.9 m Coude-Feed Telescope at Kitt Peak National 
Observatory, Arizona, USA (hereafter denoted K99) 

— 22 spectra obtained in January 2000 again using the 
2.1m telescope with the Sandiford spectrograph at 
McDonald Observatory (hereafter denoted MOO) 

Table |] shows the wavelength resolution (A/AA) at 6400 A 
and the range of signal-to-noise ratios at 6400 A for each 
dataset. The spectrograph slit widths were chosen so that 



one resolution element (FWHM of the ThAr compari- 
son lines) was 2-3 CCD-pixels. The spectra in K99 have 
considerably lower signal-to-noise ratios than the other 
datasets due to the smaller telescope. 
The spectra in S95 an d S96 were reduced using the 3A 
Software Package (Ilyin 1996 ). It uses two-dimensional dis- 
persion curves (r ows a nd columns of the CCD; see e.g. 
Duemmler et al. ( |1997 ) for a more thorough description) 
in order to calibrate the wavelength scale. The spectra in 
M99 were reduced using the 4A Software Package (Ilyin 
^000| ). The spectra in K99 and MOO were reduced us- 
ing IRAF. 4A and IRAF use three-dimensional dispersion 
curves (rows and columns of the CCD as well as time on 
the basis of several comparison spectra) to calibrate the 
wavelength scale. For all spectra, telluric lines based on 
the wavelengths given by Pierce & Breckinridge (1973) 
were used prior to heliocentric correction to establish the 
accurate wavelength zero point, correcting for tiny geo- 
metrical shifts between the comparison and stellar images 
due to bending of the spectrograph and the slit-effect, i.e. 
the shifts due to the fact that the optics is not homoge- 
neously illuminated by the stellar light (Griffin & Griffin 
197| Ilyin |2000| ) 



The projected equatorial rotational velocity vsini of the 
subgiant primary in the RS CVn binary was determined 
from our measurements using a Fourier-transform tech- 
nique described by Gray ( |l98S| ; |l992| , Ch. 17). The re- 
sult is wsini — 39kms^^. For the secondary, a value of 
vsini = 7.5kms~^ was determined by comparing spun- 
up standard spectra with that of the secondary, a value 



consistent with the one given by Vogt & Hatzes (1991). 
The radial velocities (RVs) were obtained by cross- 
correlating the UX Ari spectra with spectra of RV stan- 
dard stars. For the primary we used spectra of /3 Gem 
(KOIIIb, RV = -h(3.3±0.1)kms-^) observed in the same 
run as the UX Ari spectra, reduced in the same way and 
artificially spun up to match v sin i of the primary. For 
the secondary we used the solar FTS spectrum (Kurucz 
et al. 1984 ), artificially spun up to match v sini of the sec- 
ondary. The RVs are weighted averages of the individual 
RVs measured in several orders. The measured RVs and 
their standard deviations are given in Table 



^ Table 2 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http: / /cdsweb.u-strasbg.fr/Abstract.html . 
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In the spectra of UX Ari there are also weak lines from 
a third star present. This star was first mentioned by 
McAlister et al. (1987) when they measured the angular 
separation between the RS CVn binary and the third star. 
This and other measurements of the angular separation 
are given m Table |. From Table | it seems that the an- 
gular separation has decreased by almost 0'.'2 from 1985.8 
to 1996.9. It is not known whether the third star is part 
of the UX Ari system or a star that just happens to lie on 
the same line of sight. Its spectral classification is also not 
known, although Fabricius & Makarov (200C) determined 



its B and V magnitudes based on Hipparcos (ESA [19971) 
data. Their B and V magnitudes yield B-V = 1.19 ±0.06 
for the third star. This in turn yields a spectral type of K5 
if the star is unredden ed and on the main sequence (Gray 
1992| ). Vogt & Hatzes ( 1991 ) successfully used a synthetic 
G5 V spectrum to subtract the lines of the third star from 
the composite spectrum, and we used the solar FTS spec- 
trum as RV template to measure the radial velocities of 
this star. 

Before RV measurements of the third star could be car- 
ried out, however, the spectral flux contributions from the 
primary and the secondary had to be removed from the 
composite spectrum. Otherwise, the weak lines of the third 
star are too strongly influenced by the stronger lines of the 
other two components. 

The spectral flux contribution from the primary (sec- 
ondary) was removed using an observed spectrum of a 
single, inactive star of the same spectral classification 
as the primary (secondary). We used HD 71952 (KOIV, 
V= 6™25) and HD 84453 (KOIV, V= 6™83) for the pri- 
mary and HD 23565 (G5V, V= 7"70), HD 51419 (G5V, 
V= 6™94) and HD 71148 (G5V, V= 6™30) for the sec- 
ondary. The single star was observed in the same run as 
UX Ari and reduced in the same way. Its spectrum was 
artificially spun up to match v sin i of the primary (sec- 
ondary) , scaled to match the relative continuum flux con- 
tribution of the primary (secondary) and shifted by cross- 
correlation to the position of the primary (secondary) in 
the composite spectrum. Finally, the spun-up, scaled and 
shifted spectrum of the single star was subtracted from 
the composite spectrum. The relative continuum flux con- 
tribution for each component was determined using the 
residual line strength in the composite spectrum. 
The spectrum separation technique itself, as well as what 
applying it to our UX Ari observations can teach us about 
the three components in the UX Ari spectrum, will be de- 
scribed by Aarum & Engvold (in preparation). The results 
of applying the surface (Doppler) imaging technique to our 
UX Ari observations (and thus the details of the line pro- 
files) will be described by Aarum et al. (in preparation). 

3. The radial velocity curves of UX Ari 

3.1. The data from the literature 



A significant improvement of all orbital parameters de- 
pends strongly on the value of the period, which in turn 



is more sensitive to the overall time span covered by the 
measurements than to their quality. Thus, we supplement 
our new radial velocities with data from the literature to 
increase the time span. 

There are not many RVs of UX Ari to be found. The oldest 
dataset is due to Carlos & Popper ( |1971| ). It contains a 
few RVs obtained in the 1950s, but mostly data from 1967 
to 1970; the total number of RV pairs is 28. The second 
big dataset is due to Duquennoy et al. (1991), containing 
36 timepoints with RVs (however often only for one of 
the two stars) obtained mostly in 1977; a few RVs are 
measured in 1985-1988. Duquennoy et al. (1991) give only 
the RVs; they do not determine or improve the orbital 
parameters. Another dataset is given by Heintz (1981). 
There are only 3 RVs for each component given, and, when 
compared to a preliminary orbit, they have considerable 
scatter. This would give them such a low weight in the 
combined dataset that, together with their small number 
and the fact that their observing times overlap with the 
dataset of Duquennoy et al. (1991), we decided not to use 
them at all. 



3.2. The weights 

The optimal weights in a least squares fit are the inverse 
variances of the individual measurements. For all datasets, 
individual error estimates for the RVs are known, except 
for the set given by Carlos & Popper ( |1971| ). Thus, we 
decided to use the inverse variance as the weight, and de- 
termine an estimate for this for the Carlos & Popper set. 
An independent orbital fit of a double-lined binary RV 



curve to the data of Carlos & Popper ( 197l| ) alone was 
performed, using their relative weig htfjj. The resulting 
orbital parameters are close to those given by Carlos 
& Popper. Standard deviations for measurements having 
unit weight for the primary (1.6 km s~^) and the secondary 
(1.9kms~^) were obtained and used to calculate for each 
measurement an error by combining them with the rel- 
ative weight. These are the errors given in Table 2 for 



the measurements of Carlos & Popper (1971) and used to 



calculate the weights as the inverse variances. 
All other measurements obtained their weights as the 
inverse variances based on the errors as published. 
Preliminary orbit fits, however, indicated that a consider- 
able improvement is achieved when ignoring all measure- 
ments obtained from blended lines, i.e. with RV-differences 
< 40kms~"'^ between primary and secondary. Therefore, 
we decided to give all those measurements weight in the 
following; they are indicated by negative errors in Table 2. 
During the first orbital fits it turned out that for the 
primary is significantly larger than for the secondary, 
although the smaller accuracy of the RVs due to the much 
broader lines is already reflected by the much larger RV- 



except for the 3 very old measurements, which had weights 
in their final fit, but obtained 0.7 here; finally, since it im- 
proved the fit significantly, all data with RV-differences < 
40kms~^ between primary and secondary obtained weights 
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errors. This could be due to the systematic deviations of 
the RVs caused by the fine-structure of the hne profiles 
because of the spot activity. We therefore multiplied all 
weights of the primary by an additional factor 0.65 to 
equalize for the primary and the secondary. 

3.3. The first fit results: 7 is changing 



For ea ch da taset, i.e. Carlos & Popper ( |1971 ), Duquennoy 
et al. ( 1991 ), and each of the seasons of our own observa- 
tions, an individual orbital fit has been done. By this we 
established: 

— that there is no significant (as compared to the formal 
fit errors) change of any orbital parameter 

— except the velocity 7 of the centre of mass of the sys- 
tem, which is different for each dataset. 

The following fit, combining all datasets, thus requires 
that all orbital parameters are the same, but allows each 
dataset to obtain its own velocity zero-point, i.e. its own 
7 velocity. The results of this fit are given in Table ^ and 
compare d to the orbital parameters given by Carlos & 
Popper ( 1971 ). Here, and throughout the paper, errors 
are formal fit errors, obtained from the curvature of the 
hypersurface or from error progression, and are likely 
underestimates of the true errors. 
The results can be summarized as follows: 

— Most orbital parameters given by Carlos & Popper 
( |l97l[ ) have been confirmed. 

— Due to the much larger database and the higher ac- 
curacy of the later data, the accuracy of the orbital 
parameters has considerably increased. The period, 
in particular, is shorter than that given by Carlos & 
Popper, while still being consistent with it because of 
their larger error. 

— Despite the high significance given to the eccentric- 
ity, we do not believe that the orbit is really elliptical. 
Firstly, the eccentricity is very small; secondly, its er- 
ror is relatively large, so that the eccentricity is just a 
2.5(7 detection. Thirdly, as will be seen later, there are 
systematic effects in the radial velocity curve which 
might mimic the effects of a (small) non-zero eccen- 
tricity. Finally, the deviations from the circular orbit 
are so negligible that all other orbital parameters are 
indistinguishable from those obtained for the circular 
orbit; thus, no harm is done by neglecting e in the 
following even if it should be real. 

— The values of the velocity of the centre of mass 7 of 
the system are clearly inconsistent with each other and 
need further consideration. 

While it is not uncommon to have different velocity zero- 
points from different instruments, especially when old data 
are involved, the situation here is different. Firstly, there 
are two pairs of datasets (S95, S96 and M99, MOO) which 
were taken by the same instruments, however one year 
apart. The 7 velocities are inconsistent within the groups 



which only becomes apparent thanks to the high quality 
and number of the radial velocities leading to really small 
errors in 7. Secondly, all new data are reduced in the same 
way. This means in particular that the wavelength scale 
for each spectrum is adjusted (prior to heliocentric cor- 
rection) to the same wavelength system established by a 
large number of telluric atmospheric lines based on the 
wavelengths given by Pierce & Breckinridge (1973). This 
technique should remove all effects caused by the slit-effect 
and temporal effects like those caused by the change of am- 
bient temperature and pressure and (for the Cassegrain- 
spectrographs) bending of the spectrograph due to mo- 
tion of the telescope. Furthermore, all RVs were measured 
using the same star as template (the lAU RV-standard 
(3 Gem for the primary and the solar FTS-spectrum for 
the secondary; while (3 Gem has been observed during the 
same runs with the same instruments as UX Ari, and re- 
duced in the same way, the solar FTS-spectrum used is 
always the same). Additionally, the paper by Duquennoy 
et al. (1991) states that their radial velocities are on the 
lAU faint (my > 4.3) standard system; yet the difference 
between 7dmh and 7396 is 8.4(t. On the other hand, the 
difference between 7M99 and 7k99j which were obtained 
with different instruments a month apart, is only 1.5cr. 
The difference between 7396 and 7M00 is 44cr! Finally, it 
seems that we have, at least for the new data, a systematic 
behaviour of 7: it seems to increase from 1995 to 1996, and 
from then on it systematically decreases with time. 

For these reasons, we believe that the variation in 7 is 
genuine and not caused by any instrumental effect. The 
short-period RS CVn system UX Ari is obviously in an 
accelerated motion, most likely around the centre of mass 
with a third star. 



3.4. The final fit of the inner orbit 

If we accept this interpretation, we have to make two 
changes in order to obtain the final orbital parameters: 
First, the old datasets from the literature cover several 
years, and we see from the large difference (as compared 
to the error) between 1999 and 2000 that 7 is chang- 
ing much over one year; we thus have to subdivide all 
datasets so that each subset does not cover more than 1 
year. Secondly, if 7 is changing, so is the factor between 
the period in the observer's frame and the period in the 
restframe of the system. Given that there is no evidence 
for any change of the orbital parameters, except 7, we can 
assume that the period in the restframe of the system is 
constant, but not the one in the observer's frame. 

The final orbital fit is thus a weighted fit (using the same 
weights as before) to 17 subsets, ensuring that all parame- 
ters have the same values for all datasets, except 7. Given 
the small eccentricity found before, this fit is forced to be 
circular. The subsets have been identified in Table 2. The 
resulting orbit is given in Table ^ and shown in Fig. ^ 
Unfortunately, several of the 17 datasets consist of only 
1-3 measurements, so determination of 7 is only possible. 



R. Duemmler and V. Aarum: The radial velocities of the RS CVn star UX Ari 



5 



Table 4. Two orbital solutions using the data and errors given in Table 2; the weights of the primary have been 
additionally reduced by a factor 0.65. Fit 1: circular, 7 different 7 velocities allowed; the index of 7 identifies the 
dataset, where CP stands for Carlos & Popper ( |197lD an d DM H for Duquennoy et al. ( |199l| ). Fit 2: hke Fit 1, but 
allowing for e > 0. C&P: parameters of Carlos & Popper (1971). 



parameter 


Fit 1 


Fit 2 


C&P^ 


Pobs (days) 
Ki (kms-^) 


6.4378553 ± 0.0000046 


6.4378564 ± 0.0000046 


6.43791 ±0.00008 


57.88 ±0.17 


57.93 ±0.17 


59.4 ±0.6 


K2 (kms-i) 


66.978 ±0.033 


66.971 ±0.033 


66.7 ±0.8 


e 


0.0 (fixed) 


0.0018 ±0.0007^ 


0.0 (fixed) 


CO (deg) 




31.7 ±20.0 




To (HJD) 3 


2450642.00075 ±0.00077 


2450642.57 ±0.36 


— 


Tconj (HJD) 4 


2450640.39129 ±0.00077 


2450640.39 ±0.51 


2450640.44 ±0.13^ 


ai sin i (Rq) 


7.362 ±0.021 


7.368 ±0.021 


7.55 ±0.13^ 


a2 sin i (Rq) 
mi sin^ i {Mq) 


8.5192 ±0.0042 


8.5183 ±0.0042 


8.50 ±0.13^ 


0.6964 ±0.0066 


0.6968 ± 0.0066 


0.71 ±0.01 


7712 sin^ i (Mq) 
7CP (kms~^) 


0.6018 ±0.0055 


0.6027 ±0.0054 


0.63 ±0.01 


26.46 ±0.75 


26.47 ±0.74 


26.5 ±0.6 


7DMH (kms^-^) 


27.30 ±0.23 


27.29 ±0.23 




7S95 (kms^^) 


28.806 ±0.047 


28.911 ±0.061 




7S96 (kms~^) 


29.273 ±0.049 


29.309 ±0.053 




7M99 (kms"^) 


28.043 ±0.067 


28.094 ± 0.070 




7K99 (kms"^) 


27.898 ±0.066 


27.945 ±0.071 




7M00 (kms"-"") 


25.905 ±0.058 


25.921 ±0.061 




a (kms"^)'^ 


1.82, 0.25 


1.79, 0.25 





Note, that Carlos & Popper (1971) give m ean er rors, which have been converted to standard deviations here for consistency. 
A Lucy-Sweeney F-test (Lucy . 



Sweeney 



1971 



Lucy 1989) gives a 97.7% significance for this eccentricity. 
For e = 0, HJD of maximum RV of the primary, for e > that of periastron passage. 
HJD of the conjugation with the secondary (hotter star) in the back, 
computed from their value ("earlier star") and their period with error propagation, 
computed from their a sin i and their mass-ratio, retaining their error for a sin i. 

standard deviation of a single RV of mean weight, separately for primary and secondary, respectively. 



because no other parameter is determined for these small 
datasets alone; nevertheless, for these small datasets and 
for others due to large errors in the measurements, the 
error of 7 is sometimes so large that it hides any change 
of 7. Due to this and due to the timing of the older obser- 
vations, no indication for a changing 7 could be detected 
before our new, accurate datasets. 

In Fig. H we see that there is a nearly perfect fit of the RVs 
of the secondary. For the primary, however, we notice a 
large, systematic deviation near RV-maximum: Most new 
measurements deviate from the curve towards larger ve- 
locities. No such deviations are seen near the minimum. 
This is the reason for the larger of the primary as com- 
pared to the secondary mentioned above which motivated 
us to introduce the extra weight reduction factor of 0.65 
for the RVs of the primary. It will be very interesting to see 
the surface maps of the primary (Aarum et al., in prepa- 
ration): when the primary is moving away from us, there 
should be a strong surface feature in the facing hemisphere 
that lets us overestimate the RV, which is not visible, when 
the star shows us the other hemisphere. Also, since the 
deviation is in nearly all new measurements (which cover 
more than 4 years), this feature must be very long-lived. It 



cannot be a hot spot in the sub-secondary point, because 
this would also be visible at RV-minimum. 
This also shows that it would not be appropriate to use 
cross-correlation techniques to adjust the spectra for sur- 
face imaging: if we would have chosen a spectrum near 
RV-maximum as template and adjusted all other spectra 
to it, we would have introduced a significant wavelength 
shift. It is thus essential that the Doppler shift corrections 
prior to surface imaging come from an accurate orbit, de- 
termined over a long time span. 

3.5. The preliminary outer orbit 

Now that the 7 velocities at 17 timepoints are known, we 
could try to determine the parameters of the outer orbit. 
Unfortunately, this is not simple. While the new data show 
beyond any doubt that 7 is changing, they alone do not 
allow to determine the orbit. A maximum 7 was observed 
in 1996, but since then 7 has been decreasing steadily; no 
minimum has been observed yet. That means that from 
the new data alone, only lower limits for the period Pout 
of the outer orbit and the RV-amplitude i^out can be ob- 
tained. The old datasets, due to their unfortunate timings 
and their large errors do not improve this situation much. 
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Table 5. The final orbital solution (e = 0) using the 
data and errors given in Table 2 split into 17 datasets; 
the weights are as before. The period used is always the 
period in the restframe of the system. All parameters are 
the same for all datasets, except 7. The weighted average 
HJD for each dataset (-2400000) is given together with 
the fitted 7. 



parameter 




Prcst (days) 
Ki (kms-i) 


D.4o / z (\j6 ± u.uuuuuoy 




CV _|_ n 1 7 
/ .00 dz U.l / 




K2 (kms-i) 


DD.yoU ± U.UoD 




To (HJD) 


245Ud4z.UUzU4 zbU.UUUoi 




Tconj (HJD) ^ 


245(JD4(J.oyz7z iU.UUUoi 




ai sin i {Rq) 


7.358 ± U.U22 




02 sin i (Rq) 
mi sin^ i (Mq) 


o.oloD ± U.UU45 




0.6962 ±0.0069 




1712 sin"^ i (A^q) 
7CPi (kms"^) 


U.DUlo m U.UUOD 




27.6 ±4.0 


34785.6460 


7CP2 (kms"^) 


28.8 ±1.9 


35001.9260 


7CP3 (kms"^) 


27.61 ±0.97 


39813.4043 


7CP4 (kms~^) 


26.4 ±2.9 


39926.6130 


7CP5 (kms^^) 


28.3 ±1.1 


40125.7212 


7CP6 (kms"^) 


28.1 ±1.4 


40519.1555 


7CP7 (kms"^) 


27.5 ±1.1 


40878.2310 


7dmhi (kms"^) 


27.28 ±0.23 


43440.9205 


7DMH2 (kms^-^) 


26.86 ±0.58 


46425.6050 


7DMH3 (kms^^) 


27.67 ±0.76 


46713.0793 


7DMH4 (kms^^) 


26.77 ±0.94 


46845.1837 


7DMH5 (kms~^) 


26.62 ±0.51 


47520.6140 


7S95 (kms"^) 


28.806 ±0.049 


50055.6248 


7896 (kms~^) 


29.277 ±0.052 


50415.4416 


7M99 (kms^^) 


28.026 ±0.071 


51186.5754 


7K99 (kms^^) 


27.850 ±0.067 


51215.5638 


7M00 (kms~^) 


25.877 ±0.062 


51559.1661 


(J (kms~^)2 


1.93, 0.26 





^ HJD of the conjugation with the secondary in the back. 
^ standard deviation of a single RV of mean weight, sepa- 
rately for primary and secondary, respectively. 



ux Ari, shor: orbii {comma adjusted) 




0,0 0,2 0,'l 0,B 3 



Fig. 1. Final orbit fit (as given in Table ||) showing as 
pure errorbars the measurements from the Hterature; our 
new measurements are shown with filled circles for the sec- 
ondary (the errorbars are usually smaller than the symbol 
size) and as open circles for the primary. The different 7 
velocities have been shifted to a common value (7DMH1 ) • 



Table 6. The preliminary orbital solutions to ^U) us- 
ing the 7 values and average HJDs given in Table The 
period Pout is given in the observer's frame. 



parameter 


circular 


elliptical 


Pout (days) 
T^out (kms^^) 


3894 ± 66 


7838 ± 23 


2.90 ±0.30 


2.036 ±0.061 


e 


0.0 (fixed) 


0.622 ±0.040 


u (deg.) 




71.7±3.6 


To (HJD) 


2450495 ± 35 


2451164 ±24 


flout sin i (Rq) 


223 ±23 


247 ±13 


fim) (Mq) 
7out (kms ^) 


0.0098 ±0.0031 


0.00329 ±0.00050 


26.53 ±0.22 


27.23 ±0.12 


a (kms^^) 


0.30 


0.073 



This is with the exception of 7dhMi whose error is only 
0.23kms~^, and which will have a large impact on Pout 
and by that also allow an estimate of -ftTout- However, any 
orbital fit to the 7 velocities as they are available now has 
to be considered preliminary. 

We have done a period search in the interval Pout = 
4, 46 years. Two runs were performed, one fixing 6 = 0, 
the other allowing for an elliptical orbit. The results of 
both searches are given in Table ^; of all the orbit fits 
performed, the one yielding the smallest a is given. While 
the elliptical orbit is much better than the circular one, 
given the small number of points and their large errors we 
cannot claim it to be correct. 

The period of the circular orbit corresponds to 
(10.66 ± 0.18) yr and the projected major axis Cout sini 
to (1.04 ± 0.11) AU; the period of the elliptical or- 



bit is (21.46 ± 0.06) yr and its projected major axis is 
(1.15 ± 0.06) AU. 

3.6. The third star 

As already mentioned in Sect. ^, in the spectrum of UX Ari 
there are weak lines of a third star. Vogt & Hatzes (1991) 
measured its RV "stationary at 6.6kms~^ blueward of 
the 7 velocity of the system" . This statement refers to the 
time November 1986 to January 1987. Assuming that they 
used the 7 velocity given by Carlos & Popper ( |1971 ) , this 
velocity corresponds to 19.9 km s~^. We also measured the 
velocities of the third star in the spectrum of UX Ari after 
subtracting the spectra of the primary and the secondary. 
We did not find any variation within a run; thus, we give 
the average velocities and standard deviations in Table ^ 
It is clear that the velocity of the third star is also chang- 
ing. However, there are several reasons to believe that the 
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4. Conclusions 




Fig. 2. The preliminary orbital fits to the 17 7 velocities 
ffiven m Table |. For our new data, the last 5 points, the 
errorbars are smaller than the symbol sizes. 

Table 7. The radial velocities of the third star in the 
spectrum of UX Ari. The velocities are weighted averages 
of all RVs obtained during the corresponding run. The 
velocity from Vogt & Hatzes ( 1991 , VH86/7), referring to 
the time Nov. 1986 to Jan. 1987, is also given; the velocity 
is inferred from their description and the assumption that 



they used the 7 velocity of Carlos & Popper (1971); no 
error is available for their RV. 



dataset 


RV3 (kms-^) 


VH86/7 


19.9 


S95 


14.76±0.25 


S96 


15.532±0.086 


M99 


19.483±0.071 


K99 


19.917±0.093 


MOO 


23.43±0.42 



third star is not the body responsible for the changes of 
the 7 velocity of the short-period RS CVn system: 

— While the increase of RV3 reflects the decrease of 7, 
7 reaches its maximum in 1996, while RV3 reaches its 
minimum in (or before) 1995. 

— In 1986/7 and 1999, the third star had similar RVs. 
Either its 7 is less than the 26.5 or 27.2 kms^^ ob- 
tained from the long-period orbit of UX Ari, or its 
period is at least 24 years (it must be longer, if 19.9 
is below the 7 velocity of that binary), more than 2 
times longer than our estimate for the circular long- 
period orbit of UX Ari and still longer than the 21.5 yr 
estimated for the elliptical orbit. However, we caution 
that it is always possible that a crossing of 19.9 km s^^ 
between 1986/7 and 1995 has been missed. 

— Yet another confirmation is obtained from the conclu- 
sions presented in Sect. 0. 



By using high-accuracy radial velocities of the double- 
lined RS CVn system UX Ari we have shown that the 
7 velocity is systematically changing; all other orbital pa- 
rameters seem to be constant over time. Preliminary (cir- 
cular and elliptical) orbital solutions of 7(t) lead to a pe- 
riod of roughly 10 and 21 years, respectively. New, high- 
accuracy orbits of UX Ari are urgently needed to establish 
the accurate period of the long-period orbit. 
It is interesting that our findings do not compare well with 
those of Lestrade et al. ( 1999| ). Their angular accelerations 
measured for UX Ari, if interpreted as consequences of 
gravitational pull due to a third star, require a period of 
the outer orbit of many times their 11-year observational 
time span. Our circular orbit has a period even shorter 
than llyr; the eccentric orbit's period is less than twice 
11 yr. Furthermore, their observations obtained during JD 
w 2445000-2449000 covered, according to our eccentric 
orbit fit, the passage through the apastron, i.e. they should 
have noticed a change of direction in the proper motions. 
Our 7s are only compatible with a much larger period if 
the eccentricity is even higher than e w 0.6. Our data time 
span is not long enough to allow for such a fit. Also for 
this reason new, high-quality 7 measurements are of great 
importance. 

Several arguments (see also below) indicate that the third 
star, which produces a weak third set of lines in the spec- 
trum of UX Ari, is not the body responsible for the 7 varia- 
tions; but we found that the third star's RV is also system- 
atically varying. Thus, UX Ari is (at least) a triple system 
with a single-lined spectroscopic binary on the same line 
of sight. Also for the third star further observations are 
needed to establish its orbit and to find out whether it is 
in the fore- or in the background. 

If we take the preliminary long-period orbital solution 
seriously, the period in the observer's frame is chang- 
ing between 6'?437778 and 6'?437902 with an error of 
o-^ooooiij^ thus, the period variation is highly significant 
(11(t). Nevertheless, the phase shift caused by this varia- 
tion is negligible for Doppler imaging: the apparent motion 
of a really stationary surface feature during the 5.3 yr be- 
tween maximum and minimum is only 2°. The systematic 
change of the radial velocity by almost 6kms~"'^, however, 
is highly significant and may lead to artifacts in the maps. 
In the following, we want to draw some more conclu- 
sions about the physical parameters of the components 
in UX Ari: 

— The spectral typ e of t he secondary is given as G5V 
(Vogt & Hatzes 1991), which is consistent with the 
results from our subtraction technique. If it were of 
significantly different spectral class, the subtraction of 
G5V standards would lead to significant residuals in 
the subtracted spectrum. The classification of G5V 
implies a mass of the secondary 7712 = (0.95±0.05) Mq. 



* These numbers refer to the larger 7 variations of the circu- 
lar orbit (see Fig. B). 
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This value is obtained as an average between the val- 
ues given by Schmidt-Kaler (1982) and Gray ( 



1992 



App. B); the error accounts for the difference between 
their values as well as a small classification error. 
Together with TO2 sin i from Table | this yields an 
inclination of i = 59? 2 ± 3?3. This high inclination 
(if the orbital and rotational axes are aligned) makes 
UX Ari an ideal target for surface imaging (for results 
see Aarum et al. 1999| , Aarum et al., in preparation). 
i and mi sin'^ i (Table |^) give a mass of the primary of 
mi = (1.100 ± 0.060) Mq. 

With mi, 7712 and the mass functions /(m) from the 
outer orbit (Table |^) , and assuming that the inclina- 
tion of the long-period orbit is also 59? 2, we get for 
the mass of the third body in the system m — (0.460 ± 
0.052) Mq (for e = 0) and m = (0.307±0.020) Mq (for 
e w 0.6). The larger of these masses corresponds to an 
early M star, whose My is at least 4™ fainter than 
that of the G5 secondary, much too faint to be respon- 
sible for the third set of lines visible in the spectrum 
of UX Ari, thus confirming our earlier conclusion that 
the third star in the spectrum is a binary accidentally 
on the same line of sight. 

The same conclusion is reached in another way: If the 
third star in the spectrum would be the body responsi- 
ble for the 7 variations, the amplitude of its RV curve 
would be K3 > 26.5-14.8 = 11.7kms"^ This leads to 
a mass ratio of m^/imi + m2) = Ko^t/K^ < 0.25 and 
with the above masses mi, m2 to rris < 0.51 Mq, again 
an early M star with too small luminosity to account 
for the third set of lines. 

With the plausible assumption that the rotation of 
the two stars in the short-period system of UX Ari 
is synchronized with the orbit, i.e. Prot = -Porb = Prcst 
and irot — *orbj V sini gives for the radii of the stars 
Ri = (5.78 ± 0.13) i?0, i?2 = (1-11 ± 0.08) Rq. The ra- 
dius of the primary is consistent with its subgiant clas- 
sification; the radius of the secondary, however, is sig- 
nificantly larger than that of a G5 V star: Gray (|1992| , 
App. B) gives 0.96 i?0, while Schmidt-Kaler ( |l982|) 
gives 0.92 Rq. With such a large radius the star should 
have spectral type GO-1 V according to both authors. 
This seems to be a large deviation for the classifica- 
tion; maybe, also this star has already expanded a little 
away from the main sequence. 

With the masses and radii we obtain for the surface 
gravities loggi = 2.96 ±0.03, log 52 = 4.32 ±0.07. This 
means that the model spectra for the primary used in 
surface imaging can well be calculated with log 5 = 3; 
no interpolation between the models, given by Kurucz 
( |1993D in steps of A log 5 = 0.5, will be necessary. 
Given the inclination, we can compute a = 01+02, 
and from that and the mass ratio (Table |^) obtai n the 
effective radii of the Roche-lobes (Eggleton |l983 ), i.e. 
the radii of the spheres having the same volumes as 
the Roche-lobes. We get i?Ri,i = 7.24 ± 0.13, Rm,2 = 
6.77 ± 0.12 (the errors neglect the (at most) 1% error 
of the Eggleton ( 1983 ) approximation). 



— If we compare the Roche-lobe radii with the radii of 
the stars themselves it is clear that neither of the stars 
fills its Roche-lobe. However, the primary is close; a 
KO giant has typically a radius of lli?0 (Gray 1992) 
and since the star will go on cooling, the final radius 
should be even larger. Thus, the RS CVn system of 
UX Ari will soon become a semi-detached binary. 
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